Scattering effects on nuclear thermodynamic observables 
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The nuclear thermodynamic observables like the temperature, freeze-out volume and the specific 
heat as obtained from isotopic ratios in hot disassembled nuclear matter are examined in the light of 
the S-matrix approach to the nuclear equation of state. The values of the observables, as extracted 
without inclusion of scattering effects are found to be modified appreciably in some cases when 
scattering between the fragment species is taken care of. 
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Understanding the nuclear equation of state (EOS) is 
central in clarifying issues related to properties in com- 
pact stars, supernova dynamics or explosive nucleosyn- 
thesis P, y, 0, El ■ Laboratory experiments involving en- 
ergetic nuclear collisions are helpful to this end. The 
hot compressed nuclear system formed in these collisions 
subsequently decompresses and then becomes unstable 
against nuclear multifragmentation. It is generally be- 
lieved that this disassembly occurs at a density p consid- 
erably lower than the saturation density of normal nu- 
clear matter. It is called the freeze-out density where 
interaction among the produced fragments are assumed 
to freeze out. Relation between the different thermody- 
namic state variables like the energy E, temperature T 
and the specific volume v (=l/p) of the disassembling 
matter can be obtained from a clever manipulation of 
simplistic statistical models, broadly referred to as the 
nuclear statistical equilibrium (NSE) model that 
describe the different characteristics of the fragments pro- 
duced. 

Exploiting the NSE model Albergo, Costa, Costanzo 
and Rubbino (ACCR) Q arrived at expressions for the 
temperature and the freeze-out density in terms of the ex- 
perimental values of the double and single isotope multi- 
plicity ratios and other known quantities. The expression 
for the yield Y(A, Z) of a fragment of mass A and charge 
Z in its ground state in this model is given by 

Y{A,Z)^V^e^ A+B ^ z ^/ T , (1) 
A 

where V is the freeze-out volume, T the temperature of 
the frag menting syst em, gA the ground state degeneracy, 
\a = (\j2it / (AmT)) = A/vl, A is the nucleon thermal 
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wavelength (we use natural units, h = c = 1), [Ia is 
the chemical potential of the fragment and B(A, Z) is 
its ground state binding energy. Chemical equilibrium 
in the system assures that [ia = N/j, n + Z p, p where TV = 
A— Z and fi n and fi p are the neutron and proton chemical 
potentials. The double ratio R 2 of the fragment yields is 
then given as 

Y{A\,Z' 1 )IY{A 1 ,Z 1 ) 
Y(A' 2 ,Z 2 )/Y(A 2 ,Z 2 ) 

\AiA' 2 J 9a 1 9a' 2 

The difference in the binding energies of the selected frag- 
ments AB is given as 

AB = B(A[,Z[) - B{A U Z{) + B(A 2 , Z 2 ) - B(A' 2 , Z' 2 ).(3) 

Generally, the fragment yields Y(A, Z) are selected in 
such a way that N[ = N + 1,N^ = N 2 + 1, Z[ = Z 1} Z 2 = 
Z 2 . As T is the only unknown quantity in Eq.(2), it can 
be obtained in terms of the known quantities. 

Once the temperature T is thus known, the free neu- 
tron density p„, a measure of the volume of the fragment- 
ing system is determined from the single isotope yield 
ratio i?i = Y(A + 1, Z)/Y(A, Z) as 

Pn = Rl -7—7 ) T3 X 

\A + 1J g A +i y 
exp{[B(A,Z)-B(A + l,Z)]/T}. (4) 

The temperature obtained from Eq. (2) would yield the 
real temperature of the disassembling system provided 
the assumptions made in deriving the NSE model are 
valid. This may not be so always. One important as- 
sumption in this model is that the fragment species pro- 
duced arc taken to be noninteracting within the freeze- 
out volume. This is not true except for extremely di- 
lute systems, strong interaction corrections may modify 



2 



the fragment multiplicities and then the extracted tem- 
perature from the analysis as prescribed in Eq. (2) dif- 
fers from the real temperature of the system. Likewise, 
the free neutron density extracted through Eq. (4) us- 
ing this calculated temperature is not a true measure of 
the freeze-out volume. A number of other assumptions 
like the production of the fragments only in their ground 
states and the absence of any collective flow velocity have 
been made in deriving Eq. (1). The effects of these ap- 
proximations have already been studied earlier Q and 
would not be considered here. 

In a previous communication, it was shown that strong 
interaction corrections to the NSE model [ljl can be ap- 
propriately taken up in the S'-matrix approach [Tlj to 
the grand partition function of the dilute nuclear sys- 
tem, where in addition to all the stable mass particles, 
the two-body scattering channels between them can be 
included systematically, in principle, to all orders. This 
is the virial equation of state in terms of the S*-matrix 
elements. Assuming the dilute nuclear matter to be com- 
posed only of neutrons, protons and a-particles, Horowitz 
and Schwenk [l2j arrived at the nuclear EOS applying 
this methodology, it was extended later to include 3 H 
and 3 He [ljj. Further extension was achieved [l(| by 
including all the higher mass particles and their excited 
states as well as the scattering channels formed by any 
number of these species, estimated in the resonance ap- 
proximation. Inclusion of these heavier species is found 
to have a significant role in the evaluation of the nuclear 
EOS and other quantities [l4| like the nuclear symmetry 
energy and fragment multiplicities. Since the extraction 
of the thermodynamic state functions like the tempera- 
ture or the volume depends on the fragment multiplici- 
ties, study of the sensitivity of these extracted thermo- 
dynamic variables on the strong interaction effects in the 
nuclear medium is in order. The present communication 
aims at this study. A somewhat related study in a dif- 
ferent approach has been reported [l5| recently. 

The details of the S'-matrix (SM) approach, as applied 
to the dilute nuclear system is given in Ref. [To|. For 
completeness, a few relevant equations are presented be- 
low highlighting the approximations. All the dynamical 
information concerning the microscopic interaction in the 
system in thermodynamic equilibrium is incorporated in 
the grand partition function Z of the system, which sep- 
arates out as [11 1 



In Z = In Z, 
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The first term corresponds to contributions from all the 
nuclear species in their particle stable states behaving like 
an ideal quantum gas. The second term corresponds to 
contributions from scattering states which are absent in 
the NSE model. The first term in Eq. (5) can further be 
split into contributions from the ground and the excited 
states of the bound nuclcon clusters, so that 



The ground state contribution is given by 

d 3 p 
— r) 

Z,N " v ' 



-7 AT J V / 



(7) 



Here, the sum is over all the possible fragments includ- 
ing neutrons and protons, V is the volume of the system 
and g the ground state degeneracy of the fragments. For 
A < 8, the experimental values of the degeneracies are 
taken, for heavier nuclei they are taken to be 1 for the 
even A and 2 for odd A systems. The fragment mo- 
mentum is p and (3 is the inverse temperature 1/T. The 
effective fugacity is given by (z,n = e l3 ( fJ/Z - N+Bz - N \ Bz,n 
(= B(A, Z)) is the binding energy of the fragment and 
Hz,n{= Ha) is its chemical potential. A nucleus in a 
particular excited state is taken as a distinctly separate 
species and can be treated in the same footing as with 
the ground state. Their contributions can be obtained 
by convoluting the r.h.s. in Eq. (7) with their density 
of states in an energy interval Eq to E s where Eq is the 
first excited state taken as 2 MeV and E s is the particle 
emission threshold taken as 8 MeV. 

The scattering term in Eq. (5) can be formally written 
for nuclear matter as [Io| 

Z t .N t At & 

* e"' 3e — Tr Zt , Nt , a [AS-\e)-S(e))(8) 
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Here, the double sum refers to the sum over all possi- 
ble scattering channels, formed by taking any number of 
particles from any of the stable species (excited states in- 
cluded) and the trace is over all plane wave states for each 
of these channels. S is the scattering operator, A the bo- 
son symmctrization or the fcrmion antisymmctrization 
operator. The subscript c implies only the connected 
parts of the diagrammatics of the expression in the paren- 
theses and a denotes all other labels required to fix a 
channel in this set with proton and neutron numbers Z t 
and Nt, respectively. Bz t ,N t ,a is the sum of the individ- 
ual binding energies of all the particles in the channel 
and e is the total kinetic energy in the cm. frame of the 
scattering partners. Only two-particle scattering chan- 
nels are considered as they are expected, from binding 
energy considerations, to be more important than multi- 
particle channels with the same Z t and Nt. 

For convenience, the channels arc divided into light 
and heavy ones, consisting of low mass particles (A < 8, 
say) and heavier masses, respectively when 
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The scattering of relatively heavier nuclei is dominated 
by a multitude of resonances near the threshold; the S- 
matrix elements are then approximated by resonances 
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which, like the excited states, are again treated as ideal 
gas terms [lH [ItJ ■ Then, In Zh eavy can be written in the 
same form as \nZ ex , assuming the level density of the 
resonances to have the same functional dependence on 
energy as the excited states. We consider resonances upto 
excitations of 12 McV. The damping of the integral from 
the Boltzmann factor in Eq. (8) assures contributions 
only from relatively low energies. 

For the evaluation of In Zn g ht , only the elastic scatter- 
ing channels for the pairs NN, Nt, N 3 He, N a (N refers 
to the nucleon) and aa have been considered. These cal- 
culations involve only experimental inputs, namely, the 
phase shifts and binding energies. Details of these calcu- 
lations can be seen in pjl, [HI ■ 

Once the partition function for nuclear matter is ob- 
tained, the pressure P, the number density of the ith 
species pi, the free energy per baryon / or the entropy 
per baryon s are evaluated from the relations, 



P = T- 



MZ 
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dP 



7^ • (10) 



The total neutron and proton density are given by 



tot 
rn 



NiPi 



(11) 



with p = p l ° l + p 1 " 1 . The energy per baryon e(= / + Ts) 
and the heat capacity per baryon Cy(= (de/dT)v) can 
then be obtained. 

It is evident from the above that the fragment multi- 
plicities obtained from the NSE model are rather approx- 
imate, the presence of the scattering among the different 
fragment species modifies their multiplicities. We specif- 
ically study the influence of scattering on the extracted 
observables, namely, the temperature and the freeze-out 
volume using double and single isotope ratios and also 
on the heat capacity. For this study, we consider nuclear 
matter where Coulomb is absent and only strong inter- 
action is operative. The sums in Eqs. (7) and (8) run 
upto infinity in principle; in practice a finite sum with 
a maximum mass number A max = l000 is taken for com- 
putational facilitation. The results are found to be not 
very sensitive to further increase of A max . The binding 
energies of these nuclei are obtained using the liquid- 
drop type mass formula [l8| with Coulomb switched off. 
To obtain the effects of scattering on the aforesaid ob- 
servables, the fragment multiplicities arc calculated at a 
chosen temperature T, baryon density p and asymmetry 
x (= (Pn 1 ~ pT)Ip) in both the NSE model and the SM 
approach. One can then infer an output temperature and 
free nucleon density from the cluster composition of the 
system in the two models as prescribed in Eqs. (2) and 
(4) . One easily sees that this way in the NSE model, one 
gets back T NSE = T, the chosen temperature. The 
ratio T NSE /T SM and p™ SE / p s n M = V SM /V NSE then 



represent measure of the corrections due to scattering 
effects on the extraction of temperature and freeze-out 
volume V of the system with given p 1 T and X. The 
superscripts NSE and SM refer to quantities obtained 
in the NSE and SM models, respectively. The differ- 
ence in the inferred temperatures and energies e in the 
NSE and the SM model affects the heat capacity. Re- 
flection of this is shown in the difference between the 
specific heats at constant volume Cy M = (de/dT) v M 



and C$ SE 



(de/dT)$ SE . 

The effect of scattering on temperature from the dou- 
ble isotopic ratio is displayed in Fig. 1. Here the ratio 
of T NSE to T SM is shown as a function of the given 
temperature T at different densities, both for symmetric 
(X = 0.0) and asymmetric (X = 0.2) nuclear matter. 
The thermometer used in this calculation is the double 
isotopic ratio (t/d)/( 4 He/ 3 Hc) denoted as the t-He ther- 
mometer. The temperature T SM is always found to be 
lower than the NSE temperature; the scattering effect on 
the inferred temperature T SM increases with density. In 
general, at a given density, this effect has a peaked struc- 
ture which becomes more pronounced at lower densities. 
Isospin asymmetry is seen to be relatively more impor- 
tant only at lower temperatures, it does not affect the 
results qualitatively. So, following results are presented 
only for symmetric nuclear matter. 

Scattering effects on the extracted temperature with 
the choice of the double-ratio thermometer are presented 
in Fig. 2 at different densities as a function of tem- 
perature. Three thermometers are chosen, namely, t- 
He, (d/p)/( 4 He/ 3 He) (denoted d-He) and ( 4 He/ 3 He)/ 
( 7 Li/ 6 Li) (denoted He-Li). Choice of thermometers is 
seen to have a perceptible role on the extracted tempera- 
ture, particularly at relatively higher densities and tem- 
peratures. Unless otherwise mentioned, the subsequent 
calculations refer to t-He thermometer; among the three 
thermometers, it has the least scattering effect. 

For very dilute warm nuclear matter, it is generally 
perceived that the matter consists besides neutrons and 
protons a few light species, namely, deuterons, tritons, 
3 He and 4 He [12, LUa]- Importance of heavier clusters 
surface with increasing density and at lower temperature 
[31 . In the density and temperature range we explore, it 
is then expected that the departure of the temperature 
rpSM f rom temperature T NSE would depend on 

the selection of fragment species in the model calculation. 
In Fig. 3, we compare the ratio T NSE JT SM as a func- 
tion of temperature at different densities with the choice 
of light species upto 4 He (the light species model or LSM) 
and also with inclusion of the heavy species (the heavy 
species model or HSM). At high temperatures, the two 
results tend to merge as matter breaks up into lighter 
clusters. Significant difference between the two model 
calculations are evident at lower temperatures and rel- 
atively higher densities where heavier clusters are more 
abundant. These may be attributed to the fact that in 
the LSM, heavier species are absent by definition, conser- 
vation of baryon density increases the population of the 
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light clusters there in comparison to those in the HSM 
which introduces larger scattering effects in the LSM. 

The departure of the derived value of the freeze-out 
volume V SM from the corresponding NSE value V N 
due to strong interaction effects is shown as a function 
of temperature in Fig. 4 through the ratio y /V N 
at two relatively higher densities where the effects are 
found to be more pronounced. Model calculations have 
been performed in both HSM (upper panel) and in LSM 
(lower panel). The single isotope ratios used are (d/p) 
and (t/d) and the corresponding temperatures are ob- 
tained from the thermometers (d-Hc) and (t-He), respec- 
tively. The values of the y SM /V NSE are found to be 
sensitive to the choice of the single and double ratios. In 
LSM, the extracted volume V SM differs largely from the 
NSE value at low temperatures, but falls sharply with in- 
crease in temperature. With inclusion of heavier species, 
a qualitative change in the temperature dependence of 
ySM jyNSE is b serve H w ith a marked peak at T — 6.0 
MeV for the chosen densities. In HSM, the estimation of 
the freeze-out volume could be maximally uncertain by 
a factor of ~ 2 for the cases studied. 

Comparison is made between the values of the specific 
heat C$ M (full lines) and C* SE (dashed lines) in Fig. 5, 
where they are presented at a few densities func- 
tion of temperature inferred from the respective models. 
The calculations at p= 0.0001, 0.001, 0.01 and 0.02 far 3 
are displayed in colors red, black, magenta and cyan, re- 
spectively. The two sets of specific heats do not show any 
qualitative difference. In the S-matrix approach, because 
of scattering, the peaks in specific heat are shifted at a 
somewhat lower temperature, but their magnitudes are a 
bit higher. At a particular density, the temperature cor- 
responding to the peak in the heat capacity reflects the 
sudden onset of condensation in isochoric cooling. This is 
the condensation temperature [lj| . The scattering effects 



on the extracted thermodynamic observables are found 
to be maximum here as also can be seen from Figs. 1 and 
4(a). 

The ACCR method, in the context of the NSE model, 
gives an estimate of the thermodynamic parameters like 
the temperature and the freeze-out density of the dis- 
assembling system. They are, however, very approxi- 
mate. Various corrections are called for to improve upon 
these model values. We have here explored only the ef- 
fect of strong interactions among the different species 
produced in the nuclear disassembly on these thermo- 
dynamic observables. They involve model-independent 
entities. Significant deviations from the ACCR values are 
obtained, particularly around the condensation tempera- 
ture. The extracted temperature and the freeze-out den- 
sity are always underestimated, the scattering-corrected 
values can, however, be reconstructed either through it- 
eration or by preparing a table. 

In Ref. [lH, modifications of the ACCR values of the 
temperature and the freeze-out density are arrived at as- 
suming that the binding energies of the fragments in the 
freeze-out configuration are dressed-up due to the pres- 
ence of the medium. They considered only light frag- 
ments and their results are qualitatively the same as ours 
in the LSM. They have also incorporated effects due to 
strong interactions (they use explicit model interactions 
as opposed to ours), but the connection between the two 
approaches is not immediately obvious. 
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FIG. 1: The ratio of the temperature obtained in NSE to that 
in SM model using (t-He) thermometer at different densities 
p (in fm" 3 ) as a function of T for symmetric [X = 0.0, full 
lines) and asymmetric (X = 0.2, dashed lines) nuclear matter. 
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FIG. 2: Comparison of t /T at different densities p 
(in fm _,i ) as a function of temperature using (t-He), (d-He) 
and (He-Li) thermometers shown with full, dashed and dot- 
dashed lines, respectively. 
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FIG. 3: The same as in Fig. 1 for symmetric nuclear matter 
in HSM (full lines) and LSM (dashed lines). 
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FIG. 4: The ratio of the freeze-out volumes obtained in SM 
(V SM ) and in NSE model (V NSE ) as a function of temper- 
ature at p=0.01 and 0.02 fm -3 using the single ratios (d/p) 
and (t/d), the corresponding thermometers are (d-He) and 
(t-He). 
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FIG. 5: (color online) Heat capacities at constant volume 
calculated in SM (C£ M ) and in NSE (C$ SE ) as a function 



of T SM and T ,JVH - K , respectively at different densities. The 
red, black, magenta and cyan lines correspond to p =0.0001, 
0.001, 0.01 and 0.02 fm~ 3 , respectively. 



